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Theoretical uncertainties of (d,*He) and (*He,d) reactions due to the uncertainties of optical model 
potentials* 
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Theoretical uncertainties of single proton transfer cross sections of the (?He,d) and (d,°He) reactions 
due to the uncertainties of the entrance- and exit-channel optical model potentials are examined with the 
30Si(?He,d)*"P, 13B(d,?He)!*Be, and *4S(#He,d)?°Cl reactions at incident energies of 25 MeV, 46 MeV, and 
25 MeV, respectively within the framework of distorted wave Born approximation. Differential cross sections 
at the first peaks in the angular distributions of these reactions are found to be uncertain within around 5% due 
to the uncertainties of the optical model potentials from an result of 20000 times of calculations with the optical 
potential parameters randomly sampled. This amount of uncertainties is found to be nearly independent of the 
angular momentum transfer and the target masses within the studied range of incident energies. Uncertainties 
of the single proton spectroscopic factors obtained by matching the theoretical and experimental cross sections 


at different scattering angles are also discussed. 


Keywords: proton transfer reactions, optical model potentials, spectroscopic factors 


I. INTRODUCTION 


Proton transfer reactions, such as (d,?He) and (#He,d) reac- 
tions, are important in nuclear physics. They provide not only 
valuable nuclear reaction data for various applications, but 
also important tools for studying the single-particle structure 
of atomic nuclei, such as spectroscopic factors and asymp- 
totic normalization factors, which are of fundamental impor- 
tance in nuclear physics and nuclear astrophysics [1-8]. Re- 
action theories are necessary for extracting such structure in- 
formation from nuclear reaction measurements [9, 10]. Reli- 
able nuclear structure information relies not only on the preci- 
sion of measurements, but also on the reaction theories used. 
Experimentalists have always endeavoured to make measure- 
ments more and more precise. At the same time, it is also 
important to quantify the uncertainties of the theoretical re- 
sults of reaction cross sections [11, 12]. 

In most cases, A(d,2He)B and A(*He,d)B reactions can be 
well described with the distorted wave Born approximation 
(DWBA), in which the amplitude is expressed as: 


Taat = Sany f dra f araxS (ko, ro)ðny(R) 


nlj 


Vir f(r) XP (ka, ra), (1) 


where Yo(Ko,Ta) and xg(kg, rg) are the distorted waves 
describing the relative motions between the two particles 
in the entrance and exit channels, which are separated by 
vectors Ta and rg with wave numbers ka and kg, respec- 
tively. nı; (R) is the normalized single particle wave func- 
tion of the transferred proton in the target nucleus whose 
associated principal, orbital, and total angular momenta are 
n, land j. apj is the associated spectroscopic amplitude, 
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which is the square root of the spectroscopic factor Snij. 
f(r) = (Ys pela, r)|YalEa)) is the overlap between the in- 
ternal wave functions of ĉ°He and deuteron. The post- and 
prior-forms of the interaction Vr is Uag + Vap — Us Heg and 
Vap + Uap — Uaa, respectively, for a A(d,°He)B reaction, 
and Vap + Uaa — Uap and Uaa + Vpa — U3 Hea, respectively, 
for a A(?He,d)B reaction. Here, Upgq and Vpq are interactions 
between particles p and q, with Upq for optical model poten- 
tials, which are complex valued, and V,,, being single particle 
binding potentials, which are real. Within the framework of 
DWBA, the post- and prior-forms are equivalent. For trans- 
fer reactions with only one single particle wave function in- 
volved, the cross section is proportional to the spectroscopic 
factor, i.e., 


do \ PY 
(in) 
- do 
where (33) et 
the spectroscopic amplitude being unity. 
Experimentally, spectroscopic factors are obtained by nor- 


malizing (43) anty=1 (© the experimental data: 


do do 
P = E E l 


It is now clear that the spectroscopic factors obtained in this 
way inherit uncertainties from the uncertainties in both the 
experimental data and the theoretical calculations. 
Uncertainties in the theoretical calculations of the transfer 
reaction cross sections are rooted mainly in the uncertainties 
of the following terms: i) the optical model potentials (OMPs) 
which are responsible for the entrance and exit channel dis- 
torted waves [Xa and xg in Eq. (1)], ii) single particle poten- 
tial parameters, which are responsible for the single-particle 
wave functions [@n7; in Eq. (1)], and iii) the reaction model 
used, for instance, whether the f(r) term in the transition am- 
plitude is treated exactly or treated with zero-range approxi- 
mation, and weather the nucleon is assumed to be transferred 


do 
anlj= 


is the cross section calculated assuming 
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from one nucleus to another through a direct one-step pro- 
cess only, or through higher-order processes such as channel- 
couplings with nuclear excitations. For a given choice of re- 
action model, the uncertainties of theoretical calculations are 
mainly from uncertainties in the OMPs and the single particle 
potential parameters. 


Parameters of single particle potentials are usually deter- 
mined with the separation energy prescription, with which, 
the potential depths are determined by the separation energies 
of the transferred nucleon when the shape parameters, e.g., 
the radius parameter, ro, and diffuseness parameter, a, of a 
Woods-Saxon (WS) potential, are preselected. For WS poten- 
tials, empirical values rọ = 1.25 fm and a = 0.65 fm are fre- 
quently used [13], although attempts have been made to de- 
termine ro with, for instance, (e,e’p) measurements [14], and 
with Hartree-Fock calculations [15-18], with a fixed choice 
of a. 


Optical model potentials are usually phenomenologically 
determined by requiring them to describe elastic scattering 
angular distributions [19]. It is well-known that potential pa- 
rameters determined in this way suffer rather serious uncer- 
tainties, especially when only limited number of experimen- 
tal data were used to confine these parameters [20, 21]. This 
situation is better for systematic or global optical model po- 
tentials, whose parameters are constrained with experimental 
data, which cover rather large ranges of incident energies and 
target masses. It has been shown that usage of systematic 
OMPs helps to reduce the uncertainties of the SFs obtained 
from transfer reactions [5, 22, 23]. Due to their importance 
in the study of nuclear reactions and related subjects, such as 
nuclear structure, nuclear astrophysics, and nuclear applica- 
tions, a lot of efforts have been devoted to the study of sys- 
tematic optical model potentials [24—26]. However, most ref- 
erences of the existing systematic OMPs do not report the 
uncertainties of their parameters. Fortunately, for the OMPs 
needed in the study of (?He,d) and (d,He) reactions, the re- 
cently proposed systematic potentials of He and deuteron are 
given with uncertainties [27-29]. These uncertainties were 
obtained with the bootstrap statistical method. The boot- 
strap method simulates many repeated measurements of the 
elastic scattering data by creating new data sets of the same 
size as the original one using random sampling with replace- 
ment. Such a procedure was repeated many times, generating 
the distributions of the OMP parameters, from which the un- 
certainties of these parameters are obtained. Details of this 
method can be found in Refs. [24, 27]. These systematic 
OMPs allow us to quantify the uncertainties of (*He,d) and 
(d,*He) reaction cross sections due to the uncertainties of the 
optical model potentials. 


The reactions analyzed in this work are °°Si(?He,d)?!P 
[30, 31], '°B(d,?He)!*Be [7], and *4S@He,d)?°Cl [30, 31]. 
The angular momentum transfer, which are the same as the 
orbital angular momentum / of the corresponding single pro- 
ton wave functions with these three reactions, are l = Oh, 1h, 
and 2h, respectively. The incident energies of these reactions 
are 25 MeV, 46 MeV, and 25 MeV, respectively. These reac- 
tions are analyzed using exact finite-range DWBA taking into 
account the full complex remnant term. The single proton 
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wave function in the ground state of ĉHe is calculated with a 
p + d single particle potential provided in Ref. [32], which 
reproduces the (#He|d) overlap function calculated with the 
Green’s function Monte-Carlo method [32]. The single pro- 
ton wave functions in the ground states of target nuclei in 
the exit channels are determined with the usual separation en- 
ergy prescription with Woods-Saxon potentials whose radii 
and diffuseness parameters are rọ = 1.25 fm and ag = 0.65 
fm, respectively. The calculations are made with the com- 
puter code FRESCO [33]. 

The paper is organized as the following: the forms of op- 
tical model potentials and their parameters are introduced in 
section II, results of our theoretical calculations and the dis- 
cussion of uncertainties of the differential cross sections at 
different scattering angles are discussed in section III; the 
conclusions of this paper are summarized in section IV. 


Il. THE OPTICAL MODEL POTENTIAL PARAMETERS 


The phenomenological OMPs used in this work are defined 
as: 


U(r) = —V; fws (r)—iW, fws (r)-iW(—4au) Z fus(r)+Ve, 


(3) 
where V,, W,, and W, are respectively the depths of the real, 
volume-imaginary, and surface-imaginary parts of the central 
potential, Vc is the Coulomb potential: 


(r > Rc) 


r2 
3 RZ ; (r < Rc), 


where Zp and Zr are charge numbers of the projectile and 
targets nuclei, and Rc is the Coulomb radius of the target 
nuclei. fws is the Woods-Saxon form factor: 


ZpZre 
T ’ 


ZpZre? 
IR 


Ve(r) = (4) 


1 
~ 1 +exp(r— R;)/a;’ 


fus(r) (5) 
where R; and a; are the radius and diffuseness parameters, 
respectively with 7 = r, v, and s labeling the real term, the 
volume-imaginary term, and the surface-imaginary term in 
Eq.3. For the OMPs of *He and deuteron needed by the 
analysis of (d,?He)/(He,d) reactions, R; is calculated with 
r, Axl 3 where r; is the reduced radius parameter and Ar is 
the mass number of the target nucleus. In total, the OMPs of 
3He and deuteron projectiles both have a set of 7 OMP pa- 
rameters {V,., Tv, dv, Wy, Ws, Tw, Gw}. The uncertainties of 
these parameters can be obtained from Refs. [27, 29]. 


IHI. UNCERTAINTIES OF THE TRANSFER REACTION 


CROSS SECTIONS 


Uncertainties of the OMP parameters for He and deuteron 
for the three reactions studied in this work are listed as A, 
in Table. 1. Together shown are the uncertainties in the 
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differential cross sections caused by varying each parame- 
ter within their ranges of validity while keeping all other pa- 
rameters fixed, A,, and the uncertainties in the differential 
cross sections caused by varying all parameters simultane- 
ously, A,(total). Note that these uncertainties of OMP pa- 
rameters are given in percentage. Namely, for a parameter 
P, whose mean value being P and absolute uncertainty being 
AP, its uncertainty in percentage is: 


AP 
Ap = B x 100%. 


The absolute uncertainties of these parameters can be found 
in Ref. [27] for He and in Ref. [29] for deuteron. Un- 
certainties of the calculated differential cross sections are 
defined in the same way. The uncertainties are evaluated 
at center-of-mass angles 0em. = 0°, 9°, and 13° for the 
39Si(3He,d)*!P, '8B(d,2He)!?Be, and 34S(3He,d)’5C1 reac- 
tions, respectively, where their maximum differential cross 
sections occur. Within the ranges of these parameters, the 
cross sections depend nearly linearly on the values of each 
individual parameter. So when one parameter, p, is changed 
within the range [(1 — Ap) x p, (1+ Ap) x p] while keeping 
all other parameters fixed, the differential cross section will 
vary between its upper and lower limits, Cmax and Omin, re- 
spectively. A, (in %) is then defined as 100% x (omar — 
Omin)/ (max F Omin). 


TABLE 1. Uncertainties of the optical model potentials, Ap, their 
associated uncertainties of the differential cross sections, A,, and 
total uncertainties of the cross sections when all parameters are al- 
lowed to vary randomly, A,total, for the three reactions studied in 
this work. See the text for details. 

30Si(*He,d)*'P reaction (l = 0) 


U3He Vo Ty ay Wo Ws Tw aw 
Ap(%) 1.12 1.66 1.22 46.7 3.96 3.16 1.19 
As (%) 0.707 0.224 0.852 0.428 1.52 6.70 2.19 
Ua Vo tm Q Wy č Ws ù dw 
Ap(%) 0.964 0.151 0.129 12.2 28.6 0.390 0.134 
Ac (%) 0.669 0.036 0.064 0.377 2.39 0.405 0.057 
Ao,total (%) 4.48 
T3B(d, He) Be reaction (J = 1) 
U3He Vo Ty ay Wo Ws Tw aw 
Ap(%) 1.14 2.20 1.22 42.3 651 3.99 1.19 
Ac (%) 0.416 1.49 0.440 1.25 4.11 7.37 2.32 
Ua Vo Tv Ay Wo Ws Tw aw 
A,(%) 1.38 0.192 0.129 6.18 123 0.321 0.134 
Az (%) 0.712 0.230 0.147 1.41 0.114 0.383 0.059 
Ao, total (%) 5.27 
31S("He,d)””CI reaction (l = 2) 

U3He Vo Tu ay We Ws Tw Qw 
Ap (%) 1.11 1.56 1.22 47.7 3.97 3.15 1.19 
Az (%) 0.144 0.678 0.169 1.17 2.27 7.47 2.26 
Ua Vo Tv Ay Wo Ws Tw aw 
Ap (%) 0.960 0.151 0.129 13.7 27.1 0.400 0.134 
Az (%) 0.037 0.124 0.121 1.61 3.24 0.470 0.067 
Azv,total (%) 5.21 


Different from the fact that the differential cross sections 
depend almost linearly on each of the parameters within their 
uncertainties (when all other parameters are fixed), they tend 
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to follow a Gaussian distribution when all the OMP param- 
eters are varied simultaneously. As an example, we show in 
Fig. 1 the distributions of the differential cross sections for 
the three reactions at their first three peak angles. The cross 
sections are normalized to their mean values: 6, = (o—a)/G, 
where g are the mean values. Shown in Fig. | are results 
of 20000 times of calculations, each with the parameter sets 
{V,, Ty, Gy, Wy, Ws, Tw, Qw} randomly sampled within 
their individual uncertainties. Uniform random numbers were 
used in this procedure. This number of samples is suffi- 
cient for our analysis because the results are nearly the same 
as those with 10000 samples. Uncertainties of differential 
cross sections when all these OMP parameters simutaneously, 
A,,total, are taken to be the standard deviations of these dis- 
tributions, which are taken to be the standard deviations of 
Gaussian function that best fit these distributions in Fig. 1. 
The Ag total values for the 2°Si(?He,d)?!P, 13B(d,?He)'“Be, 
and 34S(?He,d)?°CI reactions, at scattering angles of 0°, 9° 
and 13°, respectively, are listed in Table. 1. 


From Table. 1, one observes that, 1) the uncertainties 
of OMPs in both exit and entrance channels result in an 
uncertainty of around 5% in the (*He,d) and (d,°He) reac- 
tions, which seems to be independent of the angular momen- 
tum transfer, the target masses, and the incident energies, 2) 
Aztotal ate smaller than the sums of A, for all the three re- 
actions. This suggests that some correlations exist among the 
effects of uncertainties induced by different parameters, and 
3) among all these parameters, the radius and diffuseness pa- 
rameters of the imaginary potentials, rą and aw, of the 7He 
potential, are the most sensitive to the transfer reactions. Un- 
certainties in the cross sections caused by these two param- 
eters are about twice as much as the uncertainties of these 
two parameters. This suggests that main efforts should be put 
on these two parameters in future systematic optical model 
potential studies in order to further reduce the theoretical un- 
certainties of the (*He,d) and (d,*He) reactions. 


Experimentally, spectroscopic factors studied with trans- 
fer reactions are usually obtained by matching the theoreti- 
cal cross sections to the experimental ones at the peaks of 
the angular distributions. This is a reasonable choice be- 
cause at these angles the differential cross sections are at their 
maximum values, which means that their experimental uncer- 
tainties, at least the statistical uncertainties, are the smallest. 
Spectroscopic factors obtained at these angles will have min- 
imum experimental uncertainties. However, cross sections at 
the first peaks may not always available experimentally. In 
those cases, one has to get the SFs using experimental data at 
the second or even the third peaks. It is interesting to know 
how much the uncertainties of the SFs will increase in such 
cases. This information can be find in Fig. 1 where the dis- 
tributions of cross sections at the second and the third peaks 
are given for the three reaction studied in this work, which 
have angular momentum transfer l being 0 fA, 1 A, and 2 h, 
respectively. One sees that, for the 2°Si(*He,d)?!P reaction 
which has | = Oh, the uncertainties of cross sections at the 
second and the third peaks are nearly the same, which are 
7.5% on average, larger than that with the first peak by a 
factor of 68% . On the other hand, for the '*B(d,?He)!?Be 
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Fig. 1. (Color online) Distributions of the normalized cross sections as results of calculations with 20000 times of random sampling of the 
OMP parameters at the three peaks of the angular distributions for (a) 30Si(?He,d)**P, (b) 13B(d,3He)!Be, and (c) 345(?He,d)?°Cl reactions 
at incident energies of 25, 46, and 25 MeV, respectively. The scattering angles indicated in these figures are in center-of-mass system. The 
curves represent the Gaussian functions whose heights and width parameters are found to best fit these distributions. The Y-axis is the number 
of cases whose normalized cross sections, ĝo, fall in the intervals of [6, — As, /2, ĝo + As, /2] with As, = 0.01. See the text for details. 


and 34S(?He,d)?°Cl reactions who have | = 1A and 2 A, re- 
spectively, the uncertainties of the cross sections at their first 
and second peaks are rather close, 6.0% and 5.8%, respec- 
tively, much smaller than those of their third peaks, which are 
around 11.0%. Whether these results are general for all reac- 
tions at all incident energies might be an interesting subject 
to be further studied. 


The uncertainties of theoretical cross sections at different 
scattering angles are more clearly seen in Fig. 2, in which 
the uncertainties caused by the uncertainties of the OMP pa- 
rameters are depicted as shaded bands for the three reactions 
studied in this work. The widths of these bands represent the 
upper and lower bounds of the cross sections at each scatter- 
ing angle. The corresponding uncertainties in percentage are 
shown by the brown curves with values shown by vertical or- 
dinate on the right. The following conclusions can be drown 
from this figure, 1) within their ranges of uncertainties as 
shown in Table. 1, variations of the OMP parameters mainly 
affect the amplitudes of the differential cross sections, and 
do not change the angular distributions much, especially at 
smaller angles where theoretical cross sections are compared 
with the experimental ones to get the spectroscopic factors, 
and 2) the uncertainties of differential cross sections increase 
with the increase of the scattering angles. One also observes 
that theoretical uncertainties at the shoulders of peaks may be 
even smaller than those at the peak angles. But these shoul- 
ders are where the differential cross sections change most 
abruptly with respect to the scattering angles. Spectroscopic 
factors obtained by matching the theoretical and experimental 
cross sections at these angles will have the disadvantage that 
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they will need to have much higher angular resolutions than 
those measured at the peaks. 
IV. SUMMARY 


In summary, (?He,d) and (d,*He) reactions are important 
tools to study the single particle structure of atomic nuclei. 
Knowledge of the uncertainties about theoretical calculations 
for these reactions are important for nuclear structure infor- 
mation, such as spectroscopic factors, obtained from these 
reactions. This work study the uncertainties of theoreti- 
cal cross sections for these reactions due to the uncertain- 
ties in the entrance- and exit-channel optical model poten- 
tials. Systematic potential of "He and deuteron projectiles 
are used, whose parameter uncertainties are available in Refs. 
[27] and [29], respectively. Three reactions, ?°Si(*He,d)?!P, 
13B(d,2He)!*Be, and °4S(?He,d)*°Cl, at incident energies of 
25 MeV, 46 MeV, and 25 MeV, respectively, are analyzed 
within the framework of exact finite-range DWBA. The mo- 
mentum transfer of these chosen reactions are Oh, 1A, and 
2h, respectively. The analysis is made by calculations 20000 
times for each reaction with the entrance- and exit-channel 
parameters randomly sampled simultaneously. It is found 
that uncertainties of the theoretical cross sections of these re- 
actions caused by the uncertainties of the OMP parameters 
are around 5% at scattering angles where the reactions have 
largest cross sections. Uncertainties in the single proton spec- 
troscopic factors with these reactions are concluded to be the 
same amount due to the uncertainties of the OMP parame- 
ters. Such amount of uncertainties seem to be independent of 
the angular momentum transfer and the target masses for the 
range of incident energies studied in this work. 
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2. (Color online) Uncertainty of the 30 Si He,d)’tP, 


13B(d,?He)!*Be, and 345 (?He,d)?°Cl reactions due to the uncertain- 
ties of the optical model potential parameters of both the entrance- 
and exit-channels. Widths of the bands represent upper and lower 
bounds of the cross sections (in mb/sr), which are results of 20000 
random sampling of the OMP parameters. The brown curves are 
their corresponding uncertainties in percentage, whose values are 
given by the right y-axis. 
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